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ABSTRACT 


The abilities of cloud-resolving models (CRMs) with the double-Gaussian based and the 
single-Gaussian based third-order closures (TOCs) to simulate the shallow cumuli and their 
transition to deep convective clouds are compared in this study. The single-Gaussian based TOC 
is fully prognostic (FP), while the double-Gaussian based TOC is partially prognostic (PP). The 
latter only predicts three important third-order moments while the former predicts all the third- 
order moments. A shallow cumulus case is simulated by single-column versions of the FP and PP 
TOC models. The PP TOC improves the simulation of shallow cumulus greatly over the FP TOC 
by producing more realistic cloud structures. Large differences between the FP and PP TOC 
simulations appear in the cloud layer of the second- and third-order moments, which are related 
mainly to the underestimate of the cloud height in the FP TOC simulation. Sensitivity 
experiments and analysis of probability density functions (PDFs) used in the TOCs show that both 
the turbulence-scale condensation and higher-order moments are important to realistic 
simulations of the boundary-layer shallow cumuli. 

A shallow to deep convective cloud transition case is also simulated by the 2-D versions 
of the FP and PP TOC models. Both CRMs can capture the transition from the shallow cumuli to 
deep convective clouds. The PP simulations produce more and deeper shallow cumuli than the FP 
simulations, but the FP simulations produce larger and wider convective clouds than the PP 
simulations. The temporal evolutions of cloud and precipitation are closely related to the turbulent 
transport, the cold pool and the cloud-scale circulation. The large amount of turbulent mixing 
associated with the shallow cumuli slows down the increase of the convective available potential 
energy and inhibits the early transition to deep convective clouds in the PP simulation. When the 
deep convective clouds fully develop and the precipitation is produced, the cold pools produced 
by the evaporation of the precipitation are not favorable to the formation of shallow cumuli. 
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1. Introduction 

Cloud models with higher-order turbulence closures have been used to simulate 
convective planetary boundary layer (PBL) clouds. Despite the difficulties involved in the closure 
for higher-order moments, pressure terms and dissipation terms, higher-order closure models have 
been shown to be able to simulate a wide range of boundary layers. Using a third-order closure 
(TOC) model which predicts all third-order moments and diagnoses all fourth-order moments, 
Andre et al. (1976a,b) successfully simulated penetrative convection observed in a laboratory 
experiment. Andre at al. (1978) used the same model to simulate the 24-hour evolution of the 
mean and turbulent structures of the clear convective PBL observed during the Wangara 
experiment (Clark et al. 1971). The computed daytime and nocturnal structure quantitatively 
agreed with a number of experimental data. 

Fully prognostic (FP) TOC models have been used to simulate PBL clouds with 
probability distribution function (PDF) based turbulence- scale condensation schemes (e.g. 
Sommeria and Deardorff 1977; Mellor 1977). FP TOC predicts all the third-order moments. 
Bougeault (1981a,b) investigated the effects of different turbulence-scale condensation 
parameterizations on the simulation of tradewind cumulus. A simple, skewed exponential 
distribution was found to be most efficient for parameterizing the TKE buoyancy production. 
However, the simple skewed exponential scheme does not have a variable skewness to represent 
diverse PBL clouds ranging from shallow cumulus, stratocumulus to solid stratus. Krueger (1988) 
developed a two-dimensional (2-D) cloud resolving model (CRM) that couples cloud-scale 
dynamics with a TOC that parameterizes turbulent processes in the boundary layer and in clouds. 
This TOC used a quasi-Gaussian closure (single-Gaussian based closure) to express fourth-order 
moments in terms of second-order moments and the clipping approximation to dampen the 
unrealistic growth of the triple correlations because of the use of the quasi-Gaussian closure. The 
turbulence-scale condensation scheme was based upon the single-Gaussian distribution. This 
model was used to simulate BOMEX (Barbados Oceanographic and Meteorological Experiment) 
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marine shallow cumuli (Krueger and Bergeron 1994) and was able to simulate individual shallow 
cumuli. A shortcoming of their model was seen in the simulated profile of the mean cloud 
fraction. The cloud fraction was underestimated near the cloud base and had its largest value near 
the cloud top. 

To obtain a realistic simulation of various types of PBL clouds, the internal physical 
consistency and the flexibility of TOC models must be addressed. In Bougeault’s (1981a,b) TOC 
model, the fourth-order moments were expressed by using the quasi-Gaussian closure, but 
buoyancy production terms involving liquid water were parameterized empirically based upon 
large-eddy simulation (LES) data. Recently, Lappen and Randall (2001a,b) developed a 

minimally-prognostic (MP) TOC, which predicts the triple correlation of vertical velocity ( w' ) 
only. They combined the mass flux closure (e.g. Arakawa 1969; Betts 1973; Albrecht 1979; 
Randall 1987) with a higher-order closure. A physically-based method for diagnosing the updraft 
area and the convective mass flux was developed using the second- and third-order moments of 
the vertical velocity. Golaz et al. (2002a,b; hereafter GLC) developed an MP TOC using a double - 
Gaussian based closure. There are various kinds of double Gaussians (e.g. Lewellen and Yoh 
1993; hereafter, the LY double Gaussian; Larson et al. 2002). They used the analytic double 
Gaussian I proposed by Larson et al. (2002). GLC also used a few assumptions to determine the 
magnitude of each Gaussian. All the moments in this approach are based upon the same PDL so 
that the internal physical consistency is ensured. This model was not only capable of reproducing 
reasonable profiles from first-order moments to third-order moments, but also simulating a 
reasonable time evolution of cloud fraction. Two shortcomings of the GLC model are that the 

third-order moment of liquid-water potential temperature (0^ ) was assumed to be zero and the 
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third-order moment of total water mixing ratio (qj ) was assumed to be proportional to w' . 

Cheng et al. (2004) recently compared the ability of the LP and MP TOC models to 
simulate shallow cumuli. They found that an unrealistically strong liquid water oscillation (LWO) 
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can exist in either model if the liquid-water buoyancy production terms in the third-order moment 
equations are not parameterized properly. The FP model produced a stronger LWO than the MP 
model because it has more equations with buoyancy production terms. Since the predictive 

equations of 9^' and q' t do not contain the buoyancy production term, Cheng et al. (2004) 

surmised that a partially-prognostic (PP) TOC model, which predicts w' , 0 ; ' and q} , but 

diagnoses the rest of the third-order moments, will not produce an LWO as large as that seen in 
the FP model. 

One of the purposes of this study is to evaluate the performance of a newly developed PP 
TOC model. In the PP model, w' , 0^ , and <r/ r ' are predicted and the rest of the third-order 

moments and all fourth-order moments are diagnostically obtained using a double-Gaussian PDF. 
Another purpose of this paper is to compare the abilities of CRMs with the FP and PP TOCs to 
simulate shallow cumuli and their transition to deep convective clouds. The questions that this 
paper will try to answer are 1) How do the turbulence-scale condensation and higher-order 
moments affect the shallow cumulus simulations? 2) How do the different TOCs affect the 
transition from shallow cumuli to deep convective clouds? and 3) How do the shallow cumuli and 
deep convective clouds interact? 

The rest of the paper is organized as follows. Section 2 gives a description of the PP TOC 
model. Section 3 compares the ability of the FP and PP TOCs to simulate a shallow cumulus case 
and shows the importance of the higher-order moments and the turbulence-scale condensation 
scheme with two sensitivity experiments. Section 4 explores the effects of different TOCs on the 
simulation of the transition from shallow cumuli to deep convective clouds and the interactions 
between shallow cumuli and deep convective clouds. Summary and discussions are provided in 
Section 5. 
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2. Model description 

The University of Califomia-Los Angeles (UCLA)/Langley Research Center (LaRC) 
CRM is used in this study. This model was previously known as the UCLA/Colorado State 
University CRM (Krueger 1988; Xu and Randall 1995). A PP TOC was implemented in the 
model in addition to the original FP TOC. 

The PDF is an important component of the PP version of the CRM because the higher- 
order moments are all diagnosed from the same PDF to ensure internal consistency, the cloud 
fraction and liquid water, and the buoyancy production terms of the second- and third-order 
moment equations are also derived from the PDF. Double Gaussian PDFs are able to describe the 
PBL clouds from observations (e.g., LY; Larson et al. 2002) and models (e.g., GLC). In this study, 
we used the analytical double-Gaussian II PDF proposed by Larson et al. (2002) because this PDF 
reduces to a single-Gaussian when the skewnesses of 0 / and q t go to zero. Also, unl ik e the 

analytical double-Gaussian I PDF (Larson et al. 2002), the variances of the vertical velocity of the 
two Gaussians are not required to be equal. 

In the PP model, the same first- and second-order moment equations and w' equation as 
in GLC are used except that the third-order moment equations of 0 / and q t are added in order to, 

~3 ~3 

as mentioned earlier, remove the key assumptions on 0^ and q' t used in GLC. That is, 


3 _ 3 

da' 0 ,2dw'a' 0 , ,2da T ^i ,3 dw'a' 

= 3a 3w a + XV a ^ — 

dt dz dz aaa dz 


( 1 ) 


where a can be 0/ or q { . e aaa 


10 ,3 


a' , where c 10 = 3.12 is a constant, and x is the 


9 - 1 

dissipation time scale (please see GLC for details about this time scale). X = 15 m s is the 


2 2 2 2 2 2 2 

diffusion coefficient and V“ = d~/dx + d"/dy + d"/dz ■ The terms on the right hand side are 
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the transport, the mean gradient, the dissipation, the diffusion and the fourth-order terms, 
respectively. The fourth-order term can be calculated by 


2 2 


_ 2 - - 3 


Wa' = ao wl o al (3r wal + 3w 1 a 1 + 3r wal a 1 +w l a l ) + 

2 2 ~ ~ ~ 2 ~ ~ 3 

(1 - a ) a w2°a2(3 r wa2 + 3 w 2 a 2 + 3r wa2 a 2 + w 2 a 2 ) 


( 2 ) 


where subscript 1 and 2 represent the first plume (Gaussian) and the second plume, respectively. 
a is the magnitude of the first plume, a is the standard deviation, and r is the within-plume 

correlation. Any normalized variable is calculated by Y; = (y • - x)/ o • , where x can be w , 0, 

or q , and i can be 1 or 2. Formula (2) can be derived by integrating the fourth-order moments 

over the double-Gaussian PDF. We can see that the fourth-order moments are closely related to 
the standard deviation of each plume. 


All equations are discretized on a vertically staggered C-grid. In such a grid, terms such as 
wdw' 2 /dz can be discretized without any interpolation. Turbulence advection terms, such as 
dw' 3 /dz and d w ' 4 / dz appearing in the predictive equations for w' 1 and w' 3 , respectively, can also 

be calculated directly. The production terms such as w' 2 dqJdz have to be interpolated. The 

parameters associated with the PDF are located at the mid-layer points. A detailed description of 
various elements of the model can be found in GLC, except for the aforementioned differences. 

3. Shallow cumulus case 

3. 1 Design of experiment 

BOMEX took place in June 1969. The observed cloud fraction was about 10%. It is a 
purely shallow cumulus case. The horizontal scale of most individual cumulus clouds is less than 
1 km. The cloud layer has a much lower relative humidity than that observed in stratocumulus 
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cloud layers. It is difficult for 1-D and 2-D CRMs to capture the vertical structure of shallow 
cumuli if turbulence is not properly parameterized. 

The BOMEX case is run with a vertical resolution of 40 m in a domain that extends from 
0 to 3000 m by the single column version of the UCLA/LaRC CRM. We use a time step of 2 s. 
The length of the simulation is 6 h. The surface sensible heat and latent heat fluxes are prescribed 

at 9.5 W m~ 2 and 153.4 W m~ 2 , respectively. The radiative cooling, large scale subsidence, and 
horizontal moisture advection are also prescribed. Details of the model configuration are the same 
as those described in Siebesma et al. (2003) for an LES intercomparison study. The LES- 
intercomparison data are used for comparison in this study. When the LES-intercomparison data 
are not available, those from the UCLA-LES (Stevens et al. 1999) are used. 

3.2 Results 

The computational cost of the PP model is about a half of that of the FP model, taking 66 s 
and 128 s on a 1.3 GHz notebook to run the 6 h simulations, respectively. This means that about 
one half of the computational time is used for the third-order moment equations in the FP model. 

A one-dimensional boundary-layer model needs to reasonably represent the mean profiles 
of potential temperature and water vapor mixing ratio, the distributions of cloud fraction and 
liquid water for cloudy boundary-layers, the vertical heat and moisture fluxes, the variances of 
vertical velocity, potential temperature and total water mixing ratio, the various terms in the TKE 
equation, and the skewnesses of vertical velocity, potential temperature, and total water mixing 
ratio. All these aspects are related to each other. For example, if a model underestimates the 
skewnesses of potential temperature and water vapor mixing ratio, the cloud fraction will be 
overestimated. If the TKE is underestimated, vertical transports of heat and moisture will be 
underestimated. 

The vertical distribution of clouds is an important indicator of the success of boundary 
layer models. The cloud fraction for tradewind cumulus is relatively large near the cloud base but 
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small near the cloud top (Fig. 1 a). It can be seen that the magnitude and distribution of cloud 
fraction and cloud top and base heights are quite different between the two TOC models (FigS. 
lb, c). The cloud fraction from the PP model is much closer to that of the UCLA-LES and within 
the range of LES- simulated cloud fractions (Siebesma et al. 2003). The maximum cloud fraction 
from the FP model is about 50%, which is much larger than that simulated by the UCLA-LES 
(Fig. la). Another result appearing in Fig. 1 is that the PP model simulates a quantitatively 
correct structure of tradewind cumuli, with the largest cloud fraction of about 6% near the cloud 
base, as in LESs. The FP model has its largest cloud fraction near the cloud top, which is similar 
to the result of Krueger and Bergeron (1994). This difference between the PP and FP models is 
closely related to the skewnesses of 0 / and q f , as will be discussed later. 

The differences between the PP model and the LESs are much smaller than those between 
the FP model and the LESs in all second-order moments shown in Figs. 2 and 3. The second- 

order moment of vertical velocity (wV) averaged between hours 3 and 5 of the BOMEX 
simulations are compared with the mean and standard derivations of the eleven LESs (Siebesma 
et al. 2003). There are two peaks in the LES consensus, one in the subcloud layer and the other 

within the upper part of the cloud layer, and there is a wV minimum near the cloud base between 

the two peaks. There are no significant differences of wV in the subcloud layer between the 
TOC models, but there are significant differences in the cloud layer. Although the PP model 

agrees well with the LES in the cloud layer, w'w' is underestimated significantly near the cloud 

top. The vertical profiles of 0'/ , q' 2 and q' t Qj are compared with those from the UCLA-LES 

(Figs. 2b, c, and d). Their peak magnitudes are approximately the same between the two TOC 
models. The PP model captures the vertical structure of these second-order moments better than 
the FP model because the peaks in the cloud layer are much lower in the FP model compared to 
the LES. The underestimate of the cloud layer depth in the FP model is the primary reason. The 
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vertical fluxes of liquid-water potential temperature and total water mixing ratio are shown in 
Fig. 3. It can also be seen that the differences between the PP model and the LES are much 
smaller than those between the FP model and LES. 


1 3 3 3 

The third-order moments ( w , 0^ and ql ) are better predicted in the PP model than in 

the FP model (Fig. 4). There are two w ' peaks from the LES and the two TOC models, one in 
the subcloud layer and the other in the cloud layer. The value of w ' from the FP simulation is 

about one order of magnitude smaller than the PP simulation. 0'/ and q' J from the FP model, as 

expected, are almost zero, due to the single-Gaussian based closure. The third moments produced 
by the PP model have similar magnitudes and vertical distributions to the UCLA-LES although 
some differences are apparent. This suggest that there is room for further refinement of the PP 
model. 


In summary, the PP model improves the simulations of shallow cumuli greatly over the FP 
model, in particular, reducing the magnitude of cloud fraction, correcting the unrealistic vertical 
distribution of cloud fraction and increasing the vertical extent of clouds. There are also large 
differences between the FP model and PP model for the second- and third-order moments. The 
third-order moments of 0 / , q f and w from the FP model are almost zero, while the second-order 

moments from the FP model have their peak value at lower levels than the LESs. These 
differences can be attributed to either the fourth-order moments, the turbulence- scale 
condensation or the buoyancy production terms for the higher-order moment equations since 
these are the only differences between the FP and PP models. In the following section, we will 
examine the effects of the fourth-order moments and the turbulence-scale condensation on the 
shallow cumulus simulation. The effects of the buoyancy production on the higher-order moment 
equations can be found in Cheng et al. (2004). 
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3.3 Effects of the fourth-order moments and the turbulence-scale condensation 

The FP model differs from the PP model in that it uses a single-Gaussian PDF for the 
turbulence-scale condensation, it predicts all third-order moments and diagnoses all fourth-order 
moments using the quasi-Gaussian assumption. In order to understand the importance of the PDF 
and higher-order moments in the TOC models, two experiments have been performed. In 
Experiment 1 (Exp. 1), we run the PP model except that its turbulence-scale condensation is 
replaced by the single-Gaussian formulation. In Experiment 2 (Exp. 2), we run the PP model 
while using the fourth-order moments based on the quasi-Gaussian approximation. That is, 
formula (2) is not used in this experiment, but the lower-order moments are obtained by the PP 
model, which are expected to be impacted by the inconsistent fourth-order moment closure. 

The time-height cross sections of cloud fractions from Exps. 1 and 2 are shown in Fig. 
5. Exp. 1 produces a cloud base height and a location of the maximum cloud fraction near the 
cloud base that are consistent with the UCLA-LES (Fig. la), but the cloud fraction is too small 
and the cloud top is too low. Compared with the results from the PP model (Fig. lb), it can be 
concluded that the low cloud top may be caused by the single-Gaussian condensation scheme, 
suggesting that the skewness of PDF is important near the cloud top. On the other hand, Exp. 2 
produces a cloud amount that is too high and its vertical distribution of cloud fraction is not 
similar to that of tradewind cumulus (Fig. 5b). These results are rather similar to those in the FP 
simulation (Fig. 7c). Because the common attribute between the FP simulation and Exp. 2 is the 
quasi-Gaussian based fourth-order moment closure, the poor performance of the two simulations 
seems to be related to the higher- moment closure. 

Fig. 6a shows the joint distributions of 0 / and q f diagnosed from the UCLA-LES at 780 

m between hours 4 and 5. Fig. 6b and c are the same distributions calculated from the double- 
Gaussian and single-Gaussian approaches using input from the UCLA-LES, respectively. The 
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first- and second-order moments are needed to determine the joint distribution of 0 / and q t in the 

single-Gaussian approach, while the third-order moments are also needed to determine the 
distribution in the double-Gaussian approach. This particular level is chosen because there are 
large differences in cloud fractions there among the UCFA-FES, PP model and FP model. We can 
see the effects of the PDF formulation on the turbulence-scale condensation by comparing Fig. 
6b and c. Both the single-Gaussian and the double Gaussian approaches diagnose about 10% 
cloud fractions since the PDFs are located mainly in the clear region of the diagrams (Fig. 6b, 
c). However, differences can still be seen from the diagrams. The joint distribution of 0 / and q f 

from the single-Gaussian approach is symmetric while that diagnosed from the UCFA-FES and 
that from the double-Gaussian approach are skewed. There are long tails in the cloudy region with 
joint normalized probability densities about 0.5 for both the UCFA-FES and the double-Gaussian 
approach, but greater than 1 for the single-Gaussian approach. 

Fig. 7 a and b are the same as Fig. 6 b and c, respectively, except the input is from the 
FP model. We can see that the joint distributions of 0 / and q t produced by the single-Gaussian 

and double-Gaussian approaches are generally similar. About a half of each PDF falls into the 
cloudy region of the diagrams, which means that both the single-Gaussian and double-Gaussian 
approaches predict cloud fractions of about 50%. 

Furthermore, it can be seen that the differences in cloud fraction are more pronoucedly 
caused by the higher-order moments than by the turbulence-scale condensation. This is because 
the differences between Fig. 6 and Fig. 7 are related to those in the higher-order moments, 
while those between Fig. 6 b and c, and between Fig. 7a and b are related to the turbulence- 
scale condensation. 

From the comparison above, it can be concluded that both the PDF and higher-order 
moments are important for a TOC model. Furthermore, it can be seen that the results from the FP 
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toc (Fig. 7 c) combine the unrealistic features from Exps. 1 and 2. The higher-order moments 
seem to be more important because the cloud fraction and its vertical distribution from Exp. 2 are 
more similar to those from the FP TOC than those from the PP TOC, which suggests that the 
modification of the turbulence- scale condensation scheme alone does not improve the simulation 
much. 

4. A shallow to deep convective cloud transition case 

4. 1 Design of experiment 

A shallow to deep convective cloud transition case is studied in this study. It is based on 
observations during the Large-Scale Biosphere-Atmosphere (LBA) experiment in Amazonia 
(Rondonia, Brazil) with additional measurements obtained during the TRMM-LBA (Tropical 
Rainfall Measuring Mission LBA ground validation program) on 23 February 1999. This study 
focuses on a part of the diurnal cycle from the sunrise to the early afternoon. The shallow 
convective clouds formed after a well-mixed convective boundary layer developed in the early 
morning. As the surface fluxes increased, the convective boundary layer deepened. The shallow 
cumuli then transformed to deep precipitating convective clouds around noon. 

Working Group 4 of the GCSS (GEWEX Cloud System Study; GEWEX stands for Global 
Energy and Water Cycle Experiment) chose this case to study the transition from shallow to deep 
precipitating convection (Grabowski et al., 2005). The UCLA/LaRC CRM is one of the 
participating models in this study. Grabowski et al. found that single column models (SCMs) with 
convective parameterizations produce deep convection too soon, typically less than two hours 
after sunrise. Most CRMs with 1 km horizontal grid size are capable of capturing the qualitative 
aspects of the transition from shallow to deep precipitating convective clouds, but intermodel 
differences are significant for both 2-D and 3-D simulations. 

The two-dimensional version of the UCLA/LaRC CRM with either the PP or FP TOCs is 
used to simulate this case, hereafter, the PP simulations and FP simulations. The only forcing 
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comes from the increased surface fluxes and the evolving temperature tendencies due to radiative 
processes. Both are prescribed in order to reduce the complexities caused by the surface and 
radiative feedbacks. There is no large-scale advective tendencies of temperature and moisture. 
Details of the model configuration are the same as those described in Grabowski et al. (2005). 
Briefly, the horizontal domain size is 256 km, with uniform grid size of 1 km. The vertical grid 
size is 100 m near surface and stretches to 1.5 km at 30 km. 

Two ensembles of simulations were performed. One ensemble of simulations used the PP 
TOC, while the other used the FP TOC. Each ensemble consists of five simulations, whose initial 
conditions are slightly different from each other, but whose ensemble means are the same as in 
Grabowski et al. (2005). The reason for using the ensemble of simulations is to eliminate the 
random fluctuations caused by uncertainties in the initial condition. 

4.2 Temporal evolution 

Differences between the PP and FP simulations occur when the shallow cumuli begin to 
form and become larger in the deep convective phase. The time series of the surface precipitation 
rates for the two ensembles of simulations were plotted in Fig. 8 a. Most of the precipitation is 
due to convective clouds. The precipitation produced by the PP simulations is much less than that 
produced by the FP simulations and its onset occurs about 30 minutes later. The maximum 

precipitation rate from the PP simulations is about 20 mm day" 1 and is consistent with most of the 
other CRMs participating in the intercomparison project (Grabowski et al., 2005). The shallow 
cumuli occurred near hour 2 (Fig. 8b) for both CRMs, which is consistent with all of the CRM 
experiments that used a horizontal grid size of 100 m in Grabowski et al. (2005). Clouds with 
cloud top temperatures less than 245 K occur after hour 4 (Fig. 8c). The cold cloud amount 
from the PP simulations is less than 10%, while that from the FP simulations is as large as 25%. 
The cold clouds, such as anvil or cirrus clouds, are usually associated with deep precipitating 
convective clouds. This suggests that there are more active deep convective clouds at the end of 
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the FP simulations than at the end of the PP simulations. The center of the mass of the cloud fields 


(Fig. 8d ) from the PP simulations is higher than that from the FP simulations before hour 4, 
related to deeper transitional convective clouds in the PP simulations as shown later. The total 
precipitable water increases gradually for both simulations before hour 3 (Fig. 8c) because of 
the increased positive surface total water flux. However, the total precipitable water in the FP 
simulations decreases after hour 3.5 due to the large convective precipitation, while that for the PP 
simulations still keeps increasing due to the small convective precipitation. The highest cloud 
height increases with time (Fig. 8f) when the shallow cumuli grow and some of them become 
deep convective clouds. The cloud height in the PP simulations before hour 3.5 is higher than that 
from the FP simulations, but the opposite is true when deep convective clouds develop. The 
magnitude of the difference in the cloud heights is larger in the deep convective phase than in the 
shallow cumulus phase. The liquid water path in both simulations (Fig. 8 g) increases gradually 
after hour 2, but is in a quasi-steady state at the last hour, which is consistent with the benchmark 
simulations in Grabowski et al. (2005). 

The development of deep convective clouds from shallow cumuli can be divided into 
three phases: shallow cumulus, transitional, and deep convective phases, based upon the time 
evolution of the cloud amount (Fig. 9). In the shallow cumulus phase, the cloud top is lower 
than 3 km and there are only shallow cumuli. This phase occurred before hour 3 for this case. In 
the transitional phase, the shallow cumuli grow and some of them become deep convective 
clouds. The domain-maximum cloud top height also increases from 3 km to 7 km. This phase 
occurred between hour 3 and hour 5 in the PP simulations, but between hour 3.5 and hour 4.5 in 
the FP simulations. In the deep convective phase, the deep convective clouds develop. The clouds 
can reach the tropopause. This phase occurred in the last hour of the PP simulations and lasted 
slightly longer for the FP simulations. 
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The cloud amount and its vertical structure from the PP simulations are greatly different 
from those of the FP simulations (Fig. 9). The maximum cloud amount from the FP simulations 
is comparable to the PP simulations in the shallow cumulus phase, but the cloud top is lower in 
the FP simulations. This feature was seen from the 1-D BOMEX simulation. This suggests that 
results from 1-D simulations are instructive to the more complex 2-D simulations. In the 
transitional phase, the transition from shallow cumuli to deep precipitating convective clouds 
occurs more rapidly in the FP simulations than in the PP simulations, accompanied by a rapid 
decrease of shallow cumulus population in the FP simulations. In the deep convective phase, an 
additional maximum center occurs in the upper troposphere in both simulations, but there are still 
many shallow cumuli present in the PP simulations. 

In summary, the PP simulations produce more shallow cumuli in all the phases, and 
shallow cumulus clouds are deeper in the first two phases compared with the FP simulations. In 
the deep convective phase, more precipitation and cold clouds are produced in the FP simulations 
than in the PP simulations. The onset of the precipitation in the FP simulations is about 30 
minutes earlier than in the PP simulations, due to a more rapid transition from shallow to deep 
convective clouds in the FP simulations. More detailed analyses will be presented below to 
understand the causes of these differences. 

4.3 Snapshots of cloud properties in different phases 

The cloud structure and evolution can be seen more clearly from the snapshots of cloud 

water in the three phases for representative simulations from the FP and PP ensembles (Fig. 
10 ) In both simulations, the cloud width and depth increase with time. The PP simulation 
produces more shallow clouds in all three phases, whereas a large and deep convective cloud 
occupies the right half of the domain (150 km to 200 km) for the FP simulation. It can be further 
seen that the PP simulation also produces more narrower convective clouds than the FP 
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simulation and that a spectrum of clouds with depths ranging from a couple of hundred meters to 
a couple of kilometers is present in the PP simulation. The continuous spectrum of clouds ensures 
a smooth transition from shallow cumuli to deep convective clouds because whenever a cloud 
type disappears, another cloud type grows and evolves into the type of previously dissipated 
clouds. 


The cloud-scale motion from the FP simulation is much stronger than the PP simulation. 
This can be seen from the vertical velocity fields (Fig. 11 ). In the shallow cumulus phase (not 
shown), vertical velocity is always smaller than 0.5 m s" 1 from the PP simulation, while there are 

two centers with vertical velocities greater than 0.5 m s" 1 from the FP simulation. These 
magnitudes of vertical velocities are reasonable because shallow cumuli are not associated with 
strong updrafts. In the transitional phase, vertical velocities from the PP simulation exceed 0.5 m 

s" 1 above 1 km, while the maximum centers from the FP simulation are located at much higher 

levels. There are many grid columns with vertical velocities exceeding 0.5 m s" 1 below 1 km in 
the FP simulation. This feature indicates that the re solved- sc ale transport below 1 km is more 
important in the FP simulation than in the PP simulation. In the deep convective phase, we can 
also see that vertical velocities are much larger in the FP simulation than in the PP simulation, due 
to stronger organization of convection (Fig. 1(F). 

The FP simulation produces cold pools that are stronger and occur earlier than those in the 
PP simulation. This can be seen from the perturbation field of potential temperature (Fig. 12). 
In the transitional phase, the PP simulation did not produce large potential temperature 
perturbations, whereas the FP simulation produced large areas of negative temperature 
perturbations with magnitudes less than -0.5 K below 3 km. Because the cold air inhibits the 
subgrid scale and grid-scale upward sensible heat and moisture transport, and is detrimental to the 
development of shallow cumuli, the population of shallow cumuli decreases drastically in both 
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the FP and PP simulations (Fig. 10q and f) when the cold pools are widespread in the deep 
convective phase. 

In summary, the differences in the strengths of the cold pools and cloud-scale circulations 
can explain the differences in the cloud distributions between the FP and PP simulations discussed 
earlier, especially those in the transitional and deep convective phases. 

4.4 Mean state and fluxes 

The domain- averaged potential temperature and the water vapor mixing ratio have smaller 
vertical gradients in the PP simulations than in the FP simulations (Fig. 13). They are warmer 
and drier below 1 km but colder and moister above 1 km in the PP simulations than in the FP 
simulations, respectively. Guichard et al. (2005) found that the main difference between CRMs 
and SCMs in the simulations of diurnal cycle of convection is that CRMs tend to moist and warm 
the free troposphere gradually, but the adjustments of SCMs are too quick. Fig. 13 further 
shows that the degree of moistening and warming in CRMs is highly dependent on the subgrid- 
scale parameterizations of turbulence. 

The differences in the domain-averaged states suggest that there are some differences in 
the convective instability. The time series of the convective available potential energy (CAPE) is 

shown in Fig. 14. The CAPE is calculated by CAPE = [ [g(0 c - 9 ewv )/0 ewv ]<fe, where Zq 

z o 

is the level of free convection (LFC), z 1 is the level of neutral buoyancy (LNB), g is the 

acceleration of gravity, 0 c is the potential temperature of an air parcel lifted from the surface up 

to the LFC, and 0 is the potential temperature of the environment. From Fig. 14, we can see 

that the CAPE from the FP simulations is larger than that from the PP simulations before the deep 
convective phase. The two ensembles of simulations also show a gradual increase in CAPE over 
the first 2 hours due to the increase of surface sensible and latent heat fluxes and the radiative 
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cooling. When the shallow cumuli develop, the CAPE stops increasing. It is in a quasi-steady 
state for the deep convective phase of the PP simulations, but it tends to decrease in the deep 
convective phase of the FP simulations due to the stronger stabilization effect of deep convection 
in the FP simulations. 

The differences in the domain- averaged profiles of the boundary-layer potential 
temperature and water vapor between the FP and PP simulations are closely related to the 
turbulent and resolved-scale transports. The total liquid-water potential temperature and total 
water flux (Fig. 15) from the FP simulations do not reach a level as high as that in the PP 
simulations in the shallow cumulus phase. This is consistent with the ID simulation of the 
BOMEX cumulus case, i.e., less effective turbulent transport in the FP simulations. In all the 
phases, turbulent transports in the PP simulations are much larger than those in the FP 
simulations. The turbulent transports are also much larger than the resolved- scale transports in the 
shallow cumulus and transitional phases of the PP simulations, which enhances the vertical 
mixing of the mean profiles and inhibits the increase of the CAPE. 

4.5 Interactions between shallow cumuli and deep convective clouds 

We have compared the time series, the cloud properties, the mean thermodynamic profiles 
and the resolved-scale and turbulent transports from the FP and the PP simulations. Cumulus 
clouds grow gradually and become deep convective clouds with a maximum height of 12-14 km 
in both sets of simulations. However, there are a lot of differences between the two simulations. 

The PP simulations have larger turbulent transports. The potential temperature and water 
vapor mixing ratio are more well mixed than in the FP simulations. The CAPE is smaller until 
near the end of the simulations. More and deeper shallow cumuli are produced. The deep 
convection in the PP simulations is delayed due to the stabilization effects of the shallow cumuli 
and the turbulent mixing. The productions of precipitation and cold pools are also delayed. In the 
FP simulations, the potential temperature and water vapor mixing ratio are less well mixed and 
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fewer shallow cumuli are produced. The production of the precipitation occurs earlier and the 
cold pools form earlier. The cold pools are usually located below 3 km, which is detrimental to the 
vertical transport of the sensible heat and the total water. Shallow cumuli can only be produced 
near the edge of the cold pools where the turbulent transport is not prohibited in the FP 
simulations. 

The stronger turbulent mixing associated with the shallow cumuli basically slows down 
the increase of CAPE and inhibits the early growth of the deep convective clouds in the 
transitional phases. When the deep convective clouds fully develop and precipitation is produced, 
the cold pools produced by the evaporation of the precipitation are not favorable to shallow 
cumuli, except near the edges of the cold pools. 

5. Summary and discussion 

The abilities of the double-Gaussian based and the single-Gaussian based TOC models to 
simulate the shallow cumuli and their transition to deep convective clouds are compared in this 
study. In the double-Gaussian based TOC, we use the analytical double-Gaussian II proposed by 
Larson et al. (2002). Only the third-order moments of the liquid water potential temperature, the 
total water, and the vertical velocity are predicted to determine the double-Gaussian PDF. So this 
model is called the PP model. The single-Gaussian based TOC is called the FP model since all the 
third-order moments are predicted. 

BOMEX, a “pure” shallow cumuli case, was simulated by the single-column versions of 
the FP and PP models. The PP model improves the simulation of shallow cumuli greatly over the 
FP model by producing more realistic cloud structures in both the magnitude of cloud fraction and 
the location of cloud top. Large differences between the second- and third-order moments 
produced by the FP and PP simulations appear in the cloud layer. Third-order moments of vertical 
velocity, liquid-water potential temperature and total water mixing ratio from the FP model are 
almost zero, while the second-order moments from the FP model have their peak values at lower 
levels. These differences are due to various reasons. Further experiments and analytical studies 
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have shown that both the PDF based turbulence-scale condensation and higher-order moments are 
important to realistic simulations of the boundary-layer shallow cumuli. The single-Gaussian PDF 
based turbulence-scale condensation predicts a lower cloud top, while the higher-order moments 
from the quasi-Gaussian approximation produces very large cloud faction and an unrealistic cloud 
fraction profile. 

A shallow to deep convective cloud transition case, based on observations during the LBA 
experiment in Amazonia with specific data collected during the TRMM-LBA on 23 February 
1999, was also simulated by the 2D versions of the FP and PP models. Shallow convective clouds 
transformed to deep precipitating convective clouds as the surface fluxes increased after the 
sunrise. The only forcing comes from the prescribed surface fluxes and radiative effects. There 
are no large-scale advective tendencies of temperature and moisture. These simulations have been 
analyzed in greater detail than those featured in the intercomparison study of Grabowski et al. 
(2005). 

Both the FP and PP simulations can capture qualitatively the transition from the shallow 
cumuli to deep convective clouds. However, the adjustment processes are much different. There is 
more precipitation and cold cloud produced in the FP simulations than in the PP simulations. The 
onset of the precipitation of the FP simulations is about 30 minutes earlier than that of the PP 
simulations. In the first two phases, the PP simulations produce more and deeper shallow cumuli 
than the FP simulations, but the FP simulations produce larger and wider convective clouds in the 
deep convective phase than the PP simulations. The temporal evolution of cloud and precipitation 
is found to be closely related to the turbulent transport, the cold pool and the cloud-scale 
circulation. 

The PP simulations have larger turbulent transports. The potential temperature and water 
vapor mixing ratio are more well mixed than in the FP simulations. More and deeper shallow 
cumuli are produced. The deep convection is delayed due to the stabilization effects of the 
shallow cumuli and the turbulent mixing. The production of the precipitation and the cold pool is 


20 



also delayed. In the FP simulations, the turbulent transport is much smaller. The potential 
temperature and water vapor mixing ratio are less well mixed and fewer shallow cumuli are 
produced. Both the production of the precipitation and the formation of cold pools occur earlier. 
The cold pools are usually located below 3 km, which is detrimental to the vertical transport of the 
sensible heat and the total water. Shallow cumuli can only be produced near the edges of cold 
pools in the FP simulations where the turbulent transport is not prohibited. 

The large turbulent mixing associated with the shallow cumuli slows down the increase of 
the CAPE and inhibits the early growth of the deep convective clouds in the transitional phase. 
When the deep convective clouds fully develop and the precipitation is produced, the cold pools 
produced by the evaporation of the precipitation are not favorable to shallow cumuli. 

Petch et al. (2002) studied the impact of horizontal resolution on the simulations of 
convective clouds. They reported that the precipitation, the potential temperature and total water 
fluxes are sensitive to the horizontal resolution. A sophisticated subgrid-scale scheme can reduce 
this sensitivity. We found that the differences caused by the resolution are much smaller than 
those caused by the different TOC schemes (not shown). This may be due to the fact that both 
TOC schemes can simulate the convective clouds and suggests the importance of sophisticated 
subgrid-scale schemes in reducing the sensitivity of simulations to the resolution. 

This study may shed some light on the shallow and deep convection parameterization in 
general circulation models. The turbulent transport is very important for the evolution of the deep 
convective clouds. Proper turbulent transports ensure the smooth transition from shallow cumuli 
to deep convective clouds. 
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Figure Captions 


Figure 1: Time series of (a) UCLA LES-, (b) PP-, and (c) FP-simulated cloud fraction for the 

BOMEX case. 

~~2 ~~2 

Figure 2: Profiles of (a) wV, (b) Qj~ , (c) q'~ , and (d) q' t Qj averaged between hours 3 and 5 of 

the BOMEX simulations. The solid line denotes the UCLA-LES, the dotted line denotes 
the PP model, the dashed line denotes the FP model, and the shaded area [in (a) only] 
has a width of twice the standard deviation around the mean profile calculated from 
LESs (please see Siebesma et al. 2003 for details of the LESs). 

Figure 3: Profiles of (a) p L y w'q' t , and (b) pc^vv'0'^ averaged between hours 3 and 5 of the 

BOMEX simulations, where p is the density of the air, L is the latent heat of 

evaporation and c p is the specific heat capacity. The solid line denotes the UCLA-LES, 

the dotted line denotes the PP model, the dashed line denotes the FP model, and the 
shaded area has a width of twice the standard deviation around the mean profile 
calculated from LESs. 

3 3 3 

Figure 4: Profiles of (a) w ' , (b) Qj , and (c) q' t averaged between hours 3 and 5 of the BOMEX 

simulations. The solid line denotes the UCLA-LES, the dotted line denotes the PP 
model, the dashed line denotes the FP model, and the shaded area has a width of twice 
the standard deviation around the mean profile calculated from LESs. 

Figure 5: Time series of cloud fraction for (a) EXP. 1, and (b) EXP. 2. See text for details. 

Figure 6: Joint PDFs of liquid-water potential temperature and total water mixing ratio at z = 780 

m between hours 4 and 5 of the BOMEX case produced by (a) the UCLA-LES, (b) the 
double-Gaussian PDF and (c) the single-Gaussian PDF. The moments used to 
determine the PDFs are diagnosed from the UCLA-LES. The straight line is the 
saturation curve. 
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Figure 7: 


Figure 8: 


Figure 9: 


Figure 10: 


Figure 11 : 


Figure 12: 


Figure 13: 


Figure 14: 


Joint PDFs of liquid-water potential temperature and total water mixing ratio at z = 780 
m between hours 4 and 5 of the BOMEX case produced by (a) the double-Gaussian 
PDF and (b) the single-Gaussian PDF. The moments used to determine the PDFs are 
predicted from the FP simulations. The straight line is the saturation curve. 

Time series of (a) surface precipitation, (b) total cloud fraction, (c) cold cloud fraction, 
(d) height of the center of mass of the cloud field, (e) precipitable water, (f) maximum 
cloud top height, and (g) cloud liquid water path for the shallow cumuli to deep 
convective clouds transition case. The solid line denotes the ensemble mean of the 
five PP simulations, and the dotted line denotes the FP simulations. 

Time-height cross-sections of cloud fraction for (a) PP simulations, and (b) FP 
simulations. 


Snapshots of liquid water mixing ratio (10 5 kg kg 1 ) for (a) a PP simulation at hour 
2.5, (b) an FP simulation at hour 2.5, (c) a PP simulations at hour 4, (d) an FP 
simulations at hour 4, (e) a PP simulations at hour 6, and (f) an FP simulations at hour 
6 . 


Snapshots of vertical velocity (m s' 1 ) for (a) a PP simulation at hour 4, (b) an FP 
simulation at hour 4, (c) a PP simulation at hour 6, and (d) an FP simulation at hour 6. 

Snapshots of perturbation potential temperature (K) for (a) a PP simulation at hour 4, 

(b) an FP simulation at hour 4, (c) a PP simulation at hour 6, and (d) an FP simulation at 
hour 6. 

Domain-averaged profiles of the potential temperature (a) at hour 2.5, (b) at hour 4, and 

(c) at hour 6, and the water vapor mixing ratio (d) at hour 2.5, (e) at hour 4, and (f) at 
hour 6. The solid line denotes the PP simulations, and the dotted line denotes the FP 
simulations. 

Time series of the CAPE for the FP simulations (dotted line) and the PP simulations 
(solid line). 
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Figure 15: Domain-averaged transports of the liquid water potential temperature (a) at hour 2.5, 
(b) at hour 4, and (c) at hour 6, and the total water (d) at hour 2.5, (e) at hour 4, and (f) 
at hour 6. The solid line denotes the total transport of the PP simulations, the dashed 
line denotes the turbulent transport of the PP simulations, the dotted line denotes the 
total transport of the FP simulations, and the dot-dashed line denotes the turbulent 
transport of the FP simulations. 
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